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Abstract. We consider the small time semi-classical limit for nonlin- 
ear Schrodinger equations with defocusing, smooth, nonlinearity. For a 
super-cubic nonlinearity, the limiting system is not directly hyperbolic, 
due to the presence of vacuum. To overcome this issue, we introduce new 
unknown functions, which are defined nonlinearly in terms of the wave 
function itself. This approach provides a local version of the modulated 
energy functional introduced by Y. Brenier. The system we obtain is 
hyperbolic symmetric, and the justification of WKB analysis follows. 



1. Introduction 
1.1. Presentation. We study the behavior of the solution u £ to 

(1.1) ied t u £ + —Au £ = \u £ \ 2a u £ ; uf t=0 = a £ e^ o/£ , 

as the parameter e e]0, 1] goes to zero. To fix matters, we work on M n , 
yet all the results are valid in the Torus T n . Throughout all this paper, we 
assume that the space dimension is n ^ 3, which corresponds to the physical 
cases. The unknown u £ and the initial amplitude ag are complex valued, the 
phase (fio is real- valued. The case of a more general nonlinearity, of the form 

ied t u £ + yAu e = e K f {\u £ \ 2 ) u £ ; nf t=0 = ag(x)e^/ £ , 

was discussed in [llj. In particular, WKB type analysis is justified for 
k ^ 1 (weak nonlinearity). On the other hand, when k = 0, there are only 
two cases in which the mathematical analysis of the semi-classical limit for 
nonlinear Schrodinger equations is well developed. First, for analytic initial 
data. We refer to |22] and [35J for this approach. Second, for the cubic 
defocusing nonlinear Schrodinger equation (cr = 1 in (jl.ip ). Our goal is 
to justify geometric optics in Sobolev spaces for (jl.ip when a ^ 2. This 
question has remained open since the pioneering work of E. Grenier 
where the nonlinearity had to be cubic at the origin. 
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There are several motivations to study the semi-classical limit for 
Let us mention three. First, (11. ip with cr = 2 (quintic nonlinearity) is 
sometimes used as a model for one-dimensional Bose-Einstein condensation 
(|27J). An external potential is usually considered in this framework (most 
commonly, an harmonic potential); we refer to to show that the results 
of the present paper can easily be adapted to that case (see also §6.2|) . 

Second, the limit e — ► relates classical and quantum wave equations. In 
particular, the semi-classical limit e — > for u e is expected to be described 
by the laws of hydrodynamics; see e.g. [221 1261 121|. 120] . If we assume that 
a o a o as e 0, then formally, u £ is expected to be well approximated by 
ae*^/ 6 , where 

(1.2) 



f d^ + l\V ( f ) \ 2 + \a\ 2 ' 7 



1i=0 



dt a + Vcp ■ Va + -a Ac/) = ; a\t=o = ao- 

This system is to be understood as a compressible Euler equation. Indeed, 
setting (p,v) = (\a\ 2 ,Vc/>), we see that (p,v) solves: 

( d t v + vVv + V(p (7 ) = ; v lt=0 = Vcj) , 

[ d t p + div (pv) = ; /0| t=0 = |a | 2 - 

Note that for a > 1, the above system is not directly hyperbolic symmetric, 
due to the presence of vacuum. We will see that the above system suffices to 
describe the convergences of the main two quadratic observables for u £ , that 
is, position and current densities. We will also see that passing to the limit 
e — ► in the usual conservation laws for nonlinear Schrodinger equations, 
we recover important conservation laws for the Euler equation (see §6.3|) . 
This also serves as a background to note that some blow-up results for the 
nonlinear Schrodinger equation on the one hand, and the compressible Euler 
equation on the other hand, follow from very similar identities (see §6.4j) . 
This remark reinforces the bridge noticed by D. Serre |32| . 

Another motivation lies in the study the Cauchy problem for nonlinear 
Schrodinger equations with no small parameter (e = 1 in (II. ip . typically). 
As noticed in [8, Appendix] and [9], one can prove ill-posedness results 
for energy super-critical equations by reducing the problem to semi-classical 
analysis for (II. ip . In (9j Appendix C] , a result of loss of regularity was proved 
for the cubic, defocusing nonlinear Schrodinger equation, in the spirit of the 
pioneering work of G. Lebeau [28] . It concerned the flow associated to the 
nonlinear Schrodinger equation near the origin. This was extended in |10| 
to the case of data of arbitrary size in Sobolev spaces. When the nonlin- 
earity is defocusing and not necessarily cubic, the result of [9l Appendix C] 
was extended in [2], by studying the semi-classical limit for (jl.lj) . How- 
ever, [2j does not use the complete justification of geometric optics, which 
makes it impossible to extend the results in [10] to the case of super-cubic 
nonlinearities. The analysis presented in this paper makes it possible. 
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1.2. Main results. For s ^ 0, we shall denote H s (R n ), or simply H s , 
the Sobolev space of order s. We equip H°°(R n ) = H°° := n s>0 H s (R n ) 
with the distance 



t£ 1 +\\S-9\\H' 



Note that for any k G N and any interval I, C k (I; H°°) = n s ^ C k (I; H s ). 

Assumption 1.1. We require a € N\{0} without recalling this assump- 
tion explicitly in the statements. Similarly, it is assumed that ao,4>o £ H°° , 
where recall that <f)Q is real-valued. We also suppose that Oq belongs uniformly 
to H°° and converges towards ao in H°° as e — > 0. More precisely, 

a £ = a + 0(e) in H s {R n ), Vs ^ 0. 

The first remark, based on a change of unknown due to T. Makino, 
S. Ukai and S. Kawashima |30| (see also |15| ). is that the limiting system 
(|1.2p is locally well-posed in Sobolev spaces, despite the possible presence of 
vacuum: 

Lemma 1.2 (from [2]). Let n ^ 1, and Zet Assumption li.il 6e satisfied. 
There exists T* > snc/i i/ta£ (jl.2l) /ias a unique maximal solution ((f), a) in 
C([0,T*[;F°°(]R n )). 

The proof is recalled in $2j It is based on a change of unknown intro- 
duced in |30j (see also |15j ) , which makes it possible to rewrite the equation 
under the form of a quasi-linear symmetric hyperbolic system. This trans- 
formation of the equations, which consists in introducing (v, u) := (V(/>, a a ), 
clearly exhibits a key dichotomy between a = 1 and a ^ 2. In particu- 
lar, a stability analysis in the case a ^ 2 is not straightforward, since the 
above mentioned change of variables does not seem to be well adapted to 
Schrodinger equations. 

Here is the main result of this paper. In the context of Assumption II .1\ 
we prove that the solutions of (jl.ip exist and satisfy uniform estimates on 
a time interval which is independent of e. 

Theorem 1.3. Let n ^ 3, and let Assumption \1.1\ be satisfied. There 
exists T €]0,T*[, where T* is given by Lemma such that the following 
holds. For all e G]0, 1] the Cauchy problem (jl.lfl has a unique solution 
u £ £ C ( [0 , T] ; H°° (R n )). Moreover, 

(1-3) SUp \\u £ e || Loo n T]-H k (R n )) < 

es]0,l] U ' J ' " 

where (p G C£°([0,T] x R n ) is given by (jl.2p . and the index k is as follows: 

• 1/(7 = 1, then k £ N is arbitrary. 

• If a = 2 and n = 1, then we can take k = 2. 

• If a = 2 and 2 ^ n ^ 3, then we can take k = 1. 

• If a 3 ; i/ien we can ta/«e k = a. 
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Remark 1.4. The estimate (jl.3p is trivial for = 0, from the conserva- 
tion of mass, which holds even for weak solutions ( |24j ). 

Remark 1.5. For sufficiently large a, the approach followed in this paper 
makes it possible to extend Theorem 11.31 to the case of higher dimensions, 
n ^ 4. We shall not pursue this question. 

Remark 1.6. The assumption a§ = ao + plays a crucial role in the 
above result. Indeed, the analysis in [9] shows that in the case a = 1, if we 
assume only gsq = ao + o(l), then the conclusion of Theorem 11.31 fails. For 
instance, if ag = (1 + e a )ao for some < a ^ 1, then for arbitrarily small 
t > independent of e, u £ e~ % ^l £ has oscillations of order e 1-Q \ So if a < 1, 
then u £ e~ %( ^/ £ is not bounded in if 1 . 

For a = 1, the above result is a consequence of the analysis due to 
E. Grenier |26j . and remains valid in any space dimension n ^ 1. We 
propose an alternate proof in $31 

In the quintic case a = 2, for all e > 0, the Cauchy problem (jl.ip has 
a unique global solution in C(R; H°°(W n )). Indeed, for n = 1 this follows 
from standard results for semi-linear equations; in the energy-subcritical case 
n = 2, this follows from Strichartz estimate and the conservation laws; for 
the difficult energy-critical case n = 3, this has been proved by J. Colliander, 
M. Keel, G. Staffilani, H. Takaoka and T. Tao in [19]. Therefore, the main 
point in our result is the uniform bound (jl.3|) . The same is true for the case 
n ^ 2 and a ^ 3, since the nonlinearity is then H l sub-critical. 

For a ^ 3 and n = 3, the equation is H 1 super-critical. Therefore, 
not only the bound (jl.3p is new, but also the fact that we can construct a 
smooth solution u £ to (jl.ip on a time interval [0, T] independent of e G]0, 1]. 

As a matter of fact, the proof that we present for Theorem 11.31 yields a 
stronger result, whose statement demands a more precise analysis (see The- 
orem 02]). In particular, we infer that the quadratic observables converge 
strongly towards the solution of compressible Euler equations for potential 
flows in vacuum, hence giving the Wigner measure associated to {u £ ) £ (see 
e.g. 123] for the definition and the main properties). The following result 
is proved in §4.51 

Proposition 1.7. Let n ^ 3, and let Assumption \l.l\ be satisfied. There 
exists T e]0, T*[, where T* is given by Lemma such that the position 
and current densities converge strongly on [0, T] as e — > 0: 

|n e | 2 ^|a| 2 inC([0,T];L a+1 (R n )). 

Im (eu £ Vu £ ) —> q \a\ 2 V(f) in C ([0, T];L a+1 (R n ) + L l (M n )) . 

In particular, there is only one Wigner measure associated to {u £ ) £ , and it 
is given by 

fi(t, dx, df) = \a(t, x)\ 2 dx (8) 8 (f - V<f)(t, x)) . 
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The analysis proposed to prove Theorem 11.31 allows us to compute the 
leading order behavior of the wave function u £ , provided that we know a 
more precise WKB expansion of the initial amplitude. 

Assumption 1.8. In addition to Assumption \l . 1\ we assume that there 
exists ai £ H°°(M. n ) such that 

a E = a + eoi + O (e 2 ) in H s (R n ), Vs ^ 0. 

Theorem 1.9. Let n ^ 3, and let Assumption \1.8l be satisfied. There 
exists a £ C([0, T*[; H°°), and for any T E]0,T*[, there exists e(T) > 0, 
such that u e £ C([0,T};H°°) for e €]0, e(T)], and 



(1.4) 



ae^' £ 



ae i<t>/£ 



= O(e) when a = 2 and 2 ^ n ^ 3, 

L°°([0,T];L 2 nLP) 

= 0(e) in the other cases, 

L°°([0,T];L 2 nL°°) 



where p is such that i? 1 (IR n ) C L p (M Tl ). 

Remark 1.10. In general, a ^ a, unless oq is real- valued and a\ £ iM. 
(see $5]). Therefore, the system (11.2j) does not suffice, in general, to describe 
the asymptotic behavior of the wave function u e , even though it suffices to 
describe the position and current densities (see Proposition 11.71 above). 

1.3. Scheme of the proof of Theorem 11.31 To prove that the so- 
lutions to the Cauchy problem (jl.ip exist for a time independent of e, it is 
enough to prove uniform estimates for the L°° norm of u £ (see Lemma 12.11 
below). To do so, our approach toward the semi-classical limit is to filter out 
the oscillations by the following change of unknown, involving the solution 
(a, 4>) of the limit system (II. 2|) : 

(1.5) a £ (t,x) :=u £ (t,x)e- l ^ t ' x)/£ . 

The key point is that, although it is obviously equivalent to prove L°° esti- 
mates for u £ and a £ , it is expected that one can prove uniform estimates in 
Sobolev spaces for a £ , thereby obtaining the desired L°° estimates from the 
Sobolev embedding. Obviously, uniform estimates in Sobolev spaces for u £ 
are not expected to hold, due to the rapid oscillations described by (f). 
The amplitude a £ solves the following evolution equation: 

, , ( d t a £ + V0-Va £ + -a £ A(f)-i-Aa £ = --(\a £ \ 2cT -\a\ 2a )a £ . 

(1.6) < 2 2 e V / 

[ «ft=o = a o- 

It is clear that the mass is conserved: 

\\a £ (t)\\ L2 = \\u £ (t)\\ L , = \\u £ (0)\\ L2 = \\a £ \\ L2 . 

This can be seen by multiplying (jl.6p by a £ , taking the real part and inte- 
grating over W 1 . Note that the large term in e~ 1 disappears from the energy 
estimate. Indeed, the large term in e^ 1 is a nonlinear rotation term. But 
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precisely because this term is nonlinear, it does not disappear from the esti- 
mate of the derivatives (the equation is not translation invariant). Indeed, 
Va e solves 

' d t + V</> • V + \ A<f) - 4 A) Va £ + Va e • VV<^> + \a £ VA</> 



+ - (Ja | — |a| J Va + -a V | — jot] ) = 0. 



(1.7) 

This equation is of the form 

(d t + L{v,d x ) + C{e,d x ))Va £ + i -a £ V (\a £ \ 2a - |a| 2<T ) =0, 

where C(e,d x ) is skew symmetric. Again, by multiplying (jl .T[) by Va E , 
taking the real part and integrating over R n , we obtain 

\j t W^^Wl* ~\ J div(lm(7fVa £ )) (\a £ \ 2a - \a\ 2a 

= - Re f (Va E • W(f> + \a £ V A<f) Va £ dx. 

Together with the mass conservation, this yields the following identity for 
E £ := 

ldE £ 



n 

the energy E £ := ||a e ||#i 



i / div^m^Va 2 )) (|a £ | 2<T - |a| 2<T ) ■ f,^. 



2 

for some constant depending only on the known solution (a, eft) of the 
limit system. The idea is then to find a second energy functional £ £ such 
that 

1 df e If / \ 

(1.8) 2~dt + l J div ( Im (« £Vfl£ )) (V I'" " l a | 2 J < C -A E£ + ££ )- 

By adding the two inequalities, one obtains a uniform in e energy estimate 

E £ {t) +8 £ {t) ^ e c ^\E £ (0)+£ £ (0)). 

The previous strategy has many roots. For the semi-classical limit, this goes 
back to the work of Y. Brenier [6j, P. Zhang [39 , F. Lin and P. Zhang [29J, 
and is referred to as a modulated energy estimate. Here, we will get the 
same result in a different way. Our approach amounts to trying to find a 
nonlinear change of unknown to symmetrize the equations. We will define 
g £ and q £ such that 

d t q £ + g £ div(Im(o e Va £ )) + V<£ • Vq £ + ^-^q £ A(f> = 0, 



and 



£„e _ 1 I l„£|2cr I 1 2(7 



1 9 = ~ U a I _ l a 



Not only does this allow to obtain (|1.8p with £ e := ||(? £ ||^2, but also to 
derive uniform estimates in Sobolev spaces. More precisely, we will see that 
the system of equations satisfied by (a £ ,Va £ ,q £ ) is essentially hyperbolic 
symmetric (plus some skew-symmetric terms). Therefore, we can derive 
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energy estimates, which in turn imply Theorem 11.31 Note that the idea 
of introducing new unknown functions to diminish the complexity of the 
initial problem is a strategy that has proven successful in many occasions: 
for instance, blow-up for the nonlinear wave equation, [4] (see also [3], 15]). 
low Mach number limit of the full Navier-Stokes equations [lj, or geometric 
optics for the incompressible Euler or Navier-Stokes equations [161 117|, 118] . 

2. Preliminaries 

Since, for a S N, the nonlinearity in (jl.ip is smooth, the usual theorems 
for semi-linear evolution equations (see e.g. |14| ) imply the following result. 

Lemma 2.1. Let a,n e N \ {0}. For (fixed) e e]0, 1], assume that 
uf t=0 £ H s (M. n ) with s > n/2. Then there exists T £ such that (jl.ip has 
a unique maximal solution u e £ C([0, T £ [; H s (W n )): ifT 6 < +oo, then 

(2.1) l_irn |K(i)|| Loo{Mn) = +oo. 

Consequently, ifu £ {0) G F°°(K n ), then u £ E C°°([0, T £ [; -£P°(IR n )). 

With regards to the limit system (|1.2I) . we recall the proof of Lemma [1.21 

Lemma 2.2. Let a € N andn ^ 1. For all (<p ,a ) £ H s+1 (R n )xH s (R n ) 
with s > n/2 + 1, there exists T* > such that (| 1 . 2 1) has a unique maximal 
solution (faa) in C([0, T* [; H s+1 (R n ) x fF- 1 ^)). In addition, if<p ,a £ 
H^iW 1 ), then <p,a £ C°°([0, T* [; F°°(]R n )). 

Remark 2.3. The lifespan T* is finite for all compactly support initial 
data (see Proposition EB) ■ If a = 1, then a belongs to C([0, T* [; H s (R d )) 
as soon as (</>o 5 ao) £ -£P +1 (M n ) x H s (R n ). What makes the previous result 
non-trivial is the presence of vacuum when a ^ 2: at the zeroes of a, (|1-2|) 
ceases to be hyperbolic, and this may cause a loss of regularity. 

Sketch of the proof. One can transform (|1.2j) into a quasi-linear 
system by differentiating the equation for (f>: with v = V0, one has 



(2.2) 



'd t v + v Vv + V\a\ 2cT = ; v ]t=0 = V0 O , 
c^a + w • Va + -adivu = ; a| t=0 = ao- 



2 

For the cubic case where a = 1, this system enters the standard framework 
of quasi-linear symmetric hyperbolic systems, with a constant symmetrizer. 
Thus, one can solve the Cauchy problem (|1.2p in standard fashion: one first 
solves (I2.2p and then checks that cuilv = 0, so that v = for some (p. 
By contrast, for <r > 1, System (|2.2p is no longer symmetric. However, as 
in |30| . one can prove that the Cauchy problem for (|2.2p is well-posed, with 
loss of (at most) one derivative for a, by introducing A = a a . Indeed, (v, A) 
solves a quasi- linear hyperbolic system with constant symmetrizer: 

X7\A\ 2 = ; %=o = Wo 

(2.3) 




V.4 + ^Adivu = ; A\ t=0 = ag". 
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This allows us to determine v, and hence <j>, by setting 

<t>{t,x) = Mx) ~ J (~Ht,x)\ 2 + \A(t,x)\ 2 ^J dr. 

Then dt (V(j> — v) = Vdt<j) — dtv = 0, hence v = V0. Once this is granted, 
one can define a as the solution of the second equation in (|2.2p . where v is 
now viewed as a given coefficient. Since A and a a satisfy the same linear 
equation, with identical initial data, we obtain A = a a . Therefore, (a, 0) 
solves (12. 2p . Finally, the local existence time T* may be chosen independent 
of s > n/2 + 1, thanks to tame estimates (see e.g. |34j ). □ 

For further references, we conclude this paragraph by recalling a stan- 
dard estimate in Sobolev spaces for systems of the form 

(2.4) d t U+ A j (^,U)d j U + sC(d x )U = E($,U), 

where U : [0, T) x R n -> C d with d ^ 1, e G R and: 

• $ : [0, T] x IT -> C d is a given function. 

• The Afs are dxd Hermitian matrices depending smoothly on their 
arguments. 

• C{d x ) = Ljkdjdk is a skew-symmetric second order differential 
operator with constant coefficients. 

• E a C°° function of its arguments, vanishing at the origin. 

Lemma 2.4. Let n 1 anc? s > n/2 + 1. There exists a smooth non- 
decreasing function C from [0, +oo[ to [0, +oo[ such that, for all T > 0, 
all e G R, all coefficient $ G C([0, T]; fP(M n )) and a// unknown U G 
C([0,T] ;j £P(lR n )) satisfying i/iere Zio/ds 

sup ||C/(t)||^^||C/(0)|| Hs e c ( M ) r , 
ie[0,T] 

tmtft M := ||*||_Loo([ 0) T];ifs(R™))- 

PROOF. We want to estimate the L 2 (M n ) norm of A S U, where A s is the 
Fourier multiplier (Id— A) s / 2 . To deal with smooth functions, we use the 
Friedrichs moiliners: let j G C^°(lR n ) be such that = 1 for |£| ^ 1, then 
we define J$ = j{5D x ) as the Fourier multiplier with symbol j(5£). 

With these notations, set U s := J S A S U. Since s - 1 > n/2, i? 8 " 1 ^") 
is an algebra which is stable by composition (F(u) G H s ~ 1 (M n ) whenever 
u G H'-ifR") and F G C°° satisfies F(0) =0): ?7 G C 1 ([0, T]; i2" s_2 (R n )). 
Therefore, C/ 5 is smooth: C7 5 G C 1 ([0, T]; H°°(R n )). Now write 

a t C/ 5 + A j ($,U)d j U s + eC{d x )U 6 = fs, 

with 

fs= Yl [A j (<S>,U),J s A s }d j U + J S A S E(<S>,U). 
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Since C(d x ) = -C(d x )*, and since U s G C x ([0, T]; L 2 (K n )), by taking the 
inner product in L 2 (M n ), we get 

j t \\U S \\l 2 = (d j M$,U)U s ,U s } + 2(f s ,U 5 ) 

X) ii^i(*. l7 )iii<-)ll c/ '*llL» + ii/*ii^' 

where we have used the symmetry of the matrices Aj. The Sobolev embed- 
ding and the usual nonlinear estimates (see [34] ) imply 

||5iA i (*,t0|| £oo <C(||(*,C0ll wl .<»)<^(ll(^^)llH.) J 

lp,-($,[/),j a A s ]d jt y|| i2 < if n^-^, [7)11^11^11^ <c(\\($,u)\\ HS ), 

\\J s A s E(<f>, U)\\ L2 < if COII ff . < C(||($, COII^.)i 

where Aj = Aj — Aj(0) and C denotes a smooth non-decreasing function 
independent of 5. To complete the proof, apply Gronwall lemma and let 5 
go to in the inequality thus obtained. □ 

3. Proof of Theorem 11.31 in the case er = 1 

Recall that a 6 is defined as: 

a £ {t,x) :=u £ (t,x)e- i ^ t ' x)/£ , 
where cj) £ C°°([0,T*[xR n ) is given by (jl.2p . Assume in the rest of this 



paragraph that a = 1. Then, (jl.6p reads 

d t a e + • Va £ + ^a £ A0 - i|Aa £ = -~ (ja £ | 2 - |a| 2 ) a £ . 
a ft=o = a o- 

Let s > n/2 + 1 and set t e := min(T*,T e ), where T* and T £ are given by 
Lemmas 12.11 and 12.21 We prove that there exists a function C from [0, +oo[ 
to [0, +oo[ such that, for all e g]0, 1] and all t £ [0, t £ [, 

(3.1) \\a £ (t)\\ Hs ^\\a £ (0)\\ Hs e tC ^^\ 

where 

r£ /'.A II „£| 



M e (t) :— ||a £ ||Loo([ 0)t ] ; H«(R»)) + ll( a > ( PJ\\L°°([0,t];H s + 3 (M n )) ■ 

This suffices to conclude by a standard continuity argument. Indeed, set 

M := sup ||a§[| H , + ||(a, 4>)\\l^([o,t*/2];H^<^)) < +°°> 

e6]0,l] 

and choose T g]0,T*/2] small so that M exp(T C(2M )) < 2M . Since 
M £ (0) < 2M and since M £ G C°([0,r £ [), {31} implies 

M £ (t) < 2M , Vt G [0,min{T ,T £ }[. 

Sobolev embedding then shows that ||w e (i)ll_L°c = ||a e (i)|lioo is uniformly 
bounded for all e G]0, 1] and all t G [0, min{To, T £ }[. Hence, the continuation 
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principle (pTTj) implies that T £ ^ T > for all e G]0, 1]. The estimate (fL3j) 

with <T = 1 then follows from the bound sup sup M e (t) ^ 2Mo. 

ee]0,i] t€[0,T ] 

Theorem 11.31 for a = 1 was first established by E. Grenier in [26] . whose 
approach is based on a subtle phase/amplitude representation of the solu- 
tion. Here, we give an alternate proof which consists in symmetrizing the 
large terms in e~ 1 in the equation for a £ by introducing 

I e|2 I |2 

, a — a 
q £ := 

£ 

We find directly, in view of Assumption ll.lt 

d t q £ + div(Im(a e Va £ )) + div(g e V</>) = ; \\qf t=0 \\ H °(R") = 0(1), Vs > 0. 
Furthermore, with this notation the equations for a £ and ip £ : = Va e read 

d t a £ + • Va £ + ^a e A(/> - i|Aa £ + iq £ a £ = 0, 

+ V0 • VV> £ + ^ e A(/. + V e • VV0 + ^a £ VA<p 

+ iq £ il) £ + ia £ Vq £ = i^Aijj £ . 

It is easily verified that U £ := (2q £ , a £ , a £ , ip £ ,W) G C°°([0, r e [; i7°°(]R n )) 
satisfies a system of the form ()2.4|) . that is 

+ ^ A,-($, U £ )djU £ + eL{d x )U £ = E(Q, U £ ), 

where $ := (V</>, A0, VA</>). Hence, by Lemma |2.4| we obtain the desired 
estimate (|3. 1[) and conclude the proof of Theorem 11.31 in the case a = 1. 

4. The case a ^ 2 

We now follow the strategy presented in §1.31 We introduce a nonlinear 
change of unknown functions which, together with (jl.6p . yields a quasi- 
linear system of the form (|2.4p . We conclude the proof of Theorem 11.31 and 
Proposition 1 1 . 71 thanks to a general result on the composition by non-smooth 
functions in Sobolev spaces. 

4.1. A nonlinear change of variable. As already explained, to sym- 
metrize the equations, our idea is to split the term |a e | 2<T — |a| 2<7 as a product 

\a £ \ 2 -\a\ 2 ° = g £ f3 £ = {GB)(\a £ \\\a\ 2 ) = G(r 1 ,r 2 ^ 

where (3 £ satisfies an equation of the form 

(4.1) d t p £ + L(a,4),d x )p £ + g £ div (eIm(a £ Va £ )) = 0, 

and L is a first order differential operator. Proposition 14 . 31 below shows that 
it is possible to do so. Before giving this precise statement, we introduce 
convenient notations, and explain how to formally find (3 £ . 
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Introduce the position densities 
p:= \a\ 2 £ C°°([0,T*[xR n ) ; p 6 := \a £ \ 2 = \u £ \ 2 £ C°°([0, T £ [xf n ). 
Let v = V0. Elementary computations show that: 

(4.2) <9 i/0 + div(H = 0, 

(4.3) d t p e + div Im (elfVu 6 ) = 0, 

(4.4) d t p e + div (Im(ea e Va e ) + p 6 v) = 0. 



Denote 
By writing 



J 6 :=eIm(o e Vo e ). 



ft/? 5 = (d ri B)(p £ ,p)d t p £ + (d r2 B)(p £ ,p)d t p, 

we compute, from (|4.2p and (|4.4p : 

9 t /3 £ + (5 riJ B)( / 9 £ , /3 )div(J £ + p £ W ) + (<9 r2 5)(p £ ,p) div(pu) = 0. 
Hence, in order to have an equation of the desired form (|4.1|) . we impose 

d ri B(rx,r 2 ) = G(r 1 ,r 2 ). 

Since on the other hand, 

G{r u r 2 )B(r l ,r 2 )=r!-r 2 J , 
this suggests to choose ft 6 such that 

(4.5) ((3 6 ) 2 = -^-(p 6 T +1 - 2p°p 6 + f( P ). 

cr + 1 

To obtain an operator L which is linear with respect to (3 6 we choose 

(4.6) (P 6 ) 2 = -^-(pT +1 - -^-TP a+1 - 2p%p 6 - p). 

<7 + l a + 1 

With this choice, we formally compute: 

d t l3 6 + eg 6 div(Im(cfVa £ )) + v ■ V(3 6 + divv = 0. 

Before deriving this equation rigorously, examine the right hand side of 
Taylor's formula yields 

2 (P £ T +1 - ~^-rP a+1 - 2p a (P £ ~P) = iP £ ~ P) 2 QAP £ ,P), 



a + 1 cr + 1 

where Q a is given by: 

(4.7) QAn,r 2 ) := 2a [\l - s) (r 2 + s(n - r 2 )f- 1 ds. 

Jo 

Note that there exists C a such that: 

(4.8) Q*(n,r 2 ) ^ a (rp 1 + rp 1 ) . 
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Notation 4.1. Let <reN. Introduce 



G a (n,r 2 ) = a ' = ; B a (r 1 ,r 2 ) := (n - r 2 )y/Q <7 (ri,r 2 ), 
VQ<r{ri,r 2 ) 

where Q a is given by (|4.7p and 

Example 4.2. For cr = 1,2,3, we compute 
Gi = l, B l =r l -r 2 . 

A remarkable fact is that, although the functions G a and -B CT are not 
smooth for a ^ 2, one can compute an evolution equation for the unknown 
(3 £ := £? CT (|a £ | 2 , |a| 2 ). We have the following key proposition. 

Proposition 4.3. With G a and B a as above, define 

(3 £ :=B a (\a £ \\\a\ 2 ), g £ := G a (\a £ \ 2 , |a| 2 ). 

Then(3 £ £ C x ([0, r £ [xR n ) and/ E C°([0, r £ [xR"), w/iere t £ = mm(T*,T £ ). 
Moreover, 

(4.9) d t p £ + e# £ div(Im(a £ Va £ )) + v ■ V/? £ + divw = 0. 

Remark 4.4. Again, note the dichotomy between a = 1 and a ^ 2. 
If o" = 1 then, by definition, g £ = 1 and [3 £ = p £ — p are C°° functions. 
Moreover (14. 9p simply reads 

3 t /? £ + ediv(Im(cfVa £ )) + div(v(3 £ ) = 0, 

corresponding to the equation for q £ = e^ 1 ^ 6 in Section[3l and which follows 
directly by subtracting (|4.2p from (|4.4p . 

Proof. The regularity properties of /3 £ and g £ follow from Lemmas 12.11 
and 12.21 along with the definition of f3 £ and g £ (see Notation 14. 1\ and (|4.8p ) . 

Since by definition 

p £ (d ri B a )(p £ ,p) = (p £ y- P °, 

P £ (d r2 B a )(p £ ,p) = o-(p°-p°- 1 p £ ), 
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we have 

(i £ d t (5 e 

= (3 £ (d ri B a )(p £ ,p)d tP £ + (3 £ (d r2 B (T )(p £ ,p)d t p 

= -p £ (d ri B a )(p £ ,p) div(J e + p £ v) - f3 e (d r2 B a )(p e ,p) div(pv) 
= -((PT - P a ) div( J e + P £ v) - a{p° - p°~ x p £ ) div(H- 
From this we compute 

fi £ (d t f3 £ + eg £ div(Im(a £ Va £ )) + v ■ V(3 £ + ^^P £ div = 0. 



Introduce 



uj £ :={p £ =p} = { (t, x) £ [0, T £ [xR n | p £ (t, x) = p(t, x) } 
= ([0,r £ [xM n )\{/3 £ /0} (by (48])). 



Then gj]) holds on ([0,r e [xM n ) \ lo £ , and hence on ([0,r £ [xR")\w £ by 
continuity. To prove the proposition, it thus suffices to show 

d t p £ + eg £ div(Im(o e Va £ )) + v ■ V/3 £ + ^y^/? £ div?; = on u £ , 

where A denotes the interior of the set A. Since (3 £ = on uj £ , it is enough 
to prove that div(Im(a £ Va £ )) = on uj £ . This in turn follows from (|4.2p 
and (|4.4p . which yield: 

&\v{lm{lfVa £ )) = -e' 1 (d t (p £ - p) + div((p £ - . 

This completes the proof. □ 



We will see that \a £ \ 2a — \a\ 2a is of order 0(e), so we naturally set 
^ £ := Va £ ; q £ := e~ 1 p £ . 
We infer from the previous computations that (a £ ,ip £ ,q £ ) solves: 



(4.10) < 



-ig £ q £ a £ . 



dta £ + v ■ Va e + ^a £ divv — *^Aa £ 

d t ^ £ + v ■ X7ip £ + -iff div v + ip £ ■ Vv + t^ £ V div u - i^Aip 6 
= -iq £ V (a £ g £ ) - ia £ g £ Vq £ , 

d t q £ + v ■ Vq £ + g £ div(lm(a £ -0 e )) + —^—q £ divu = 0. 



Simply by writing 



g £ div(Im(a £ ^ £ )) = lm(g £ af div^> £ 
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we can rewrite the previous system as 



(4.11) 



d t a £ + v • Va £ 



i-Aa £ 
2 



1 



-a £ divv — ig £ q £ a e 



d t ij) e + v ■ Vip £ + ia £ g £ Vq £ - i-Aij) £ = 

1 1 

ip £ div v — tp £ • Vt> — -a e V div v — iq £ V (a £ g £ 



2 T T 2 

d t q £ + v ■ Vq £ + lm(g £ a £ divip £ ) 

Note that in view of Assumption ll.il 



a + l 



g £ divv. 



(4.12) 



l |t=o 



fP(R n ) 



+ 



^1 



t=0 



_fP(R n ) 



O(l), Vs > 0. 



A similar estimate for the initial data of q £ is a more delicate issue, since 
B a is not a smooth function. We postpone this estimate to §4.21 

The left hand side of (14. IIP is a first order quasi-linear symmetric hyper- 
bolic system, plus a second order skew-symmetric term. The right hand side 
can be viewed as a semi-linear source term. We deduce from Proposition l4.3l 



Corollary 4.5. On [0,r E [xIR n , the function U £ := (2q £ , a £ , a £ , i> £ , tp ) 
satisfies an equation of the form 

8 t U £ + A ] {v,a £ g £ ,lfg £ )d ] U £ + eL{d x )U £ = E{^,U £ ,a £ g £ ,V{a £ g £ )), 

where $ = (V0, V 2 </>, V 3 0), the Aj's are Hermitian matrices linear in their 
arguments, C{d x ) = ^Lj^djdk is a skew- symmetric second order differ- 
ential operator with constant coefficients, and E is a C°° function of its 
arguments, vanishing at the origin. 

Theorem 11.31 and Proposition 11.71 are consequences of the following: 

Theorem 4.6. Let n ^ 3, and let Assumption li.il be satisfied. There 
exists T e]0,T*[, where T* is given by Lemma such that the following 
holds. For all e G]0, 1], the Cauchy problem (II. ip has a unique solution 
u £ E C([0,r];#°°(IR n )). Moreover, 



(4.13) 



sup 

ee]o,i] 



\L°°([0,T];H k (R n )) 



+ \\<f 



\L°°([0,T];H k 



< +00, 



where the index k is as follows: 

• If a = 1, then k S N is arbitrary. 

• If a = 2 and n = 1, then we can take k = 2. 

• If a = 2 and 2 ^ n ^ 3, then we can take k = 1. 

• If a ^ 2>, then we can take k = a. 

Considering the first part of the estimate (|4.13j) yields Theorem 11.31 We 
explain why Proposition 11.71 is a consequence of (|4.13p in §4.51 
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4.2. Quasi- linear analysis. We now have to estimate (a £ ,q £ ,ip £ ) in 

Sobolev spaces. Let us briefly explain the difficulty. To clarify matters, 
suppose that g £ = G(|a e | 2 , |a| 2 ) for some smooth function G G C°°(IR 2 ). In 
particular this is so in the cubic case a = 1. Then, in view of Corollary 14.5^ 
U £ := (2q £ ,a £ ,a £ ,Tp £ ,iF) G C([0, r e [; H°°(R n )) 3+2n satisfies a system of the 
form U23D, 

d t U £ + Aj(^,U £ )djU £ + eL(d x )U £ = £($,U £ ), 

where := (|a| 2 , Vp, V0, V 2 0, V 3 0). The key difference with the system 
in Corollary 14.51 is the absence of dependence upon the extra unknown g £ . 
Then, Lemma 12.41 yields estimates in Sobolev spaces (of arbitrary order). 

Assume now a ^ 2. One can check that the previous symmetrization 
provides us with uniform a priori estimates in l? . However, the estimates 
of the derivatives require a careful analysis. Indeed, recall that 

g £ = G a (\a £ \ 2 ,\a\ 2 ) with G ff (n,r 2 ) - P " ( '' L -'" 2) 



\/Qcr (n,r 2 ) 

where P a and Q a are defined in Notation 14.11 and (|4.7p respectively. There- 
fore, G a need not be smooth at the origin. The classical approach, which 
consists in differentiating the equations, thus certainly fails here. Yet, as we 
will see, we need only estimate a £ g £ in H° . Introduce 

(4.14) F (T {z,z') = zG (T (\z\ 2 Az'f) : a £ g £ = F a (a £ , a) . 



One can check that F a G C CT_1 but F a G" C a . Hence, to estimate a £ g £ in 
H a , one cannot use the usual nonlinear estimates. Instead, we will use that 
F a is homogeneous of degree a and the following lemma. 

Lemma 4.7. Let p ^ 1 and m ^ 2 be integers and consider F: MP — > C. 
Assume that F G C°°(IR P \ {0}) is homogeneous of degree m, that is: 

F(Xy) = X m F(y), VA > 0, Vy G M p . 

Then, for n ^ 3, there exists K > such that, for all u G ff m (IR n ) wii/i 
vaZueu in W, F(u) G F m (M n ) and 

H-F'(lt) Hjym ^ K \\u\\ff m ■ 

The same is true when m = 1 and n G N. 

Remark 4.8. Note that the result is false for n ^ 4 and m ^ 2. Also, 
one must not expect -F(m) G -ff m+1 (M™), even for u G iT°°(R n ). For instance, 
if 

.2 

n = 1 = p, m = 2, -F(y) = y |y| , ■u(x) = xe x , 

then -F(u) G H 2 (R) and -F(u) i/ 3 (IR). Similarly, in general, one must not 
expect F a (u,v) G H a+1 (R n ), even for (u,v) G ^(M™) 2 . 
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PROOF. We prove the result by induction on m. Consider first the case 
m = 2. Observe that, by assumption, F G C m ~ 1 (W). To regularize F, let 

X € Cfi° (W) be such that < x < 1> x(l/) = 1 for \v\ < 1 and x(j/) = for 
\y\ ^ 3, with |Vx(y)| < 1. For £ G N, define F e G C°°(W) by 

F,(y) = (l- X (^))F(y). 

We claim that, for all y £W and all £ G N, 

|i^(y)| < C F |y| 2 , < 4C F \y\ , |0A*K*/)l < 4C F , 

where dj = d yj and 

Cp := sup \F(z)\ + sup sup |<9jF(z)| + sup sup \djd k F(z)\ . 

|z|<3 ' l<j<p|z|<3 l<j,fc<p|z|=l 

Since vanishes in a neighborhood of the origin, it suffices to establish 
these bounds for y ^ 0. The first bound follows from the homogeneity: 
\Fi{y)\ ^ \F(y)\ = \y\ 2 \F(y /\y\)\. For the second one, compute 

djF e (y) = (1 - X ^y))9jF(y) - r\d jX ){£y)F{£y), 

where we used £F(y) = £~ l F(£y). Since 1 ^ \£y\ ^ 3 on the support of 
(djx)(£y), and since djF: MP — > C is homogeneous of degree 1, we infer 

\djF e (y)\ ^\y\ ( sup \djF (z)\ + 3 sup \(d jX )(z)F(z)\ ) ^4C F \y\. 

\\z\^3 2GRP J 

The same raisoning yields 

djd k F e (y) = (1 - X {iy))djdkF{y) - (d jX )(£y)d k F(£y) 
- {d k x){£y)d 3 F(£y) - (djdkxWyMiy). 

The last three terms are clearly bounded by Cp since \£y\ ^ 3 on the support 
of x(ty)- Also, the first term is bounded by Cp since djd k F: MP \ {0} — > C 
is homogeneous of degree 0. This completes the proof of the claim. 

With these preliminary established, we easily obtain that there exists K 
such that for all £ G N and all u G H 2 (M n ) with values in W p , 

\\F e (u)\\ L2 ^K\\u\\ Loo \\u\\ L2 , 

||Vfi(u)|| La ^K\\u\\ Loo \\Vu\\ L2 , 

||V 2 F^(u)|| L2 ^ K\\u\\ LOO ||V 2 u|| L2 + i^T ||Vw||^ 4 . 

The Sobolev embeddings H^M 11 ) C L 6 (M n ) and H 2 (R n ) C L°°(IR n ) for 
n G {1,2,3} then imply that there exists a constant K such that, for all 
£ G N and all u G H 2 (M n ), 

\\F e (u)\\ H 2^K\\u\\ 2 H2 . 
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This in turn implies the desired result for F(u) by using the dominated 
convergence theorem and a duality argument. Indeed, for all <p £ Cq°(W 1 ), 



F(u)<pdx 



lim / Fi(u)<pdx 
^+00 J 



< lim sup \\F e (u)\\ H2 |M|#-2 



< K lkllfi-2 IMI#-2 » 

which implies F{u) £ H 2 (W l ) together with 11^(^)11^2 «S K \\u\\ 2 H2 . 

Assume now the result at order m ^ 2, and prove the result at order 
m+ 1. Let F £ C°°(M. P \ {0}) be homogeneous of degree m+ 1. We have 

||-F(U)|| L 2 ^ K ||«||™oo |M| L 2 < • 

Since m > 3/2 n/2, i2" m (M n ) is an algebra and 

\\VF(u)\\ Hm ^ K \\Vu\\ Hm \\F'(u)\\ Hm . 

By assumption, F' £ C°°(IR P \ {0}) is homogeneous of degree m, hence the 
induction assumption yields: 

ll-F'Ml.ffm ^ k \\ u \\h™ ■ 

Therefore, 

\\VF(u)\\ Hm ^KWuW^tli. 
The case m = 1 can be treated in a similar fashion. □ 

The lemma turns out to be useful to estimate the source term in (|4.11|) . 
but also to estimate the initial data for q £ . By definition, we have 

1 2 I-/I2 



(z,z')=(a e ,o)' 



where 



Q a (z, z') = yjQ a {\z\*. | "-' 



'12^ 

1/2 



2a 



J\l-s){\z'\ 2 + s{\z\ 2 -\z'\ 2 )y- 1 ds\ 



The function Q a is not smooth, but homogeneous of degree a — 1. So when 
(7 3, we can estimate g £ in H a ~ l at time i = thanks to this lemma. See 

S3 

To complete the proof of Theorem 14.61 i n view of Lemma 12.14 we seek 
an H 2 estimate of a 6 , since 

H 2 (R n ) C L°°(M. n ), n < 3. 

This boils down to an H 1 estimate of U e defined in Corollary 14.51 How- 
ever, the estimates of the derivatives require a careful analysis. Indeed, the 
classical approach, which consists in differentiating the equations, certainly 
fails here because G a need not be smooth. Moreover, Lemma 12.41 requires 
to control U e in H s with s > n/2 + 1, so we would demand s = 3 for n = 3 
and s £ N. In view of Lemma 12.41 and Corollary 14.51 we have to estimate 
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a £ g £ in H 4 . Because of the lack of smoothness of G a , such an estimate 
seems hopeless in general. We therefore proceed in two steps. First, using 
the particular structure exhibited in Corollary 14.51 we relax the assumption 
s > n/2 + 1 in Lemma f2.41 to s > n/2. Next, we use Lemma f4.7l to overcome 
the lack of smoothness of G a , and obtain the desired a priori estimates. 

Proposition 4.9. Assume a ^ 2. Let U £ be the vector-valued function 
given by Corollary \4-5\ and m > n/2. Then for all for t £ [0,t £ [, it satisfies 
the following a priori estimate: 

sup \\U%s)\\ Hm ^ \\U £ (0)\\ Hm e tc ^ Me ^\ 
se[CM] 

With M £ (t) := \\^\\L°°{[Q y t};H m ) + ll^ £ |L°°([0,t];if m ) + H a£ S' £ |lL°°([0,t];H m + 1 )- 

Sketch of the proof. Resume the proof of Lemma |2T4"1 The quan- 
tities that appear in M £ are those on the last three lines of the proof of 
Lemma 12.41 First, we have: 

\\VA j (v,a £ g £ ,a £ g £ )\\ Loo ^ C (\\v\\ wl>00 + \\a £ g E \\ w i, x ) 

^ C (IMI#m+i + \\a £ g £ \\ Hm+1 ) . 

Next, since Aj is linear in its arguments, Aj = Aj and: 

II I A? ' A m ]<9j :U £ ||^2 Si K \\Aj \\fj m + l ||f^ £ ||jym 

< C(\\v\\ Hm+1 + \\a £ g £ \\ )\\u £ \\ 

Finally, 

\\E(^,U £ ,a £ g £ ,V(a £ g £ ))\\ Hm ^ C (\M Hm ,\\U £ \\ Hm ,\\a £ g £ \\ Hm+1 ) . 
We conclude the proof thanks to Gronwall lemma. □ 

4.3. The case a > 3. Recall that from (gj3|), 

a £ g £ = F a (a £ ,a) , 

where F a is homogeneous of degree a. For a ^ 3 and n ^ 3, Lemma 14.71 
yields 

\\a £ g e \\ H ° ^K(\\a £ \\ H . + \\a\\ H ,) a . 

Proposition 14.91 with m = a — 1 ^ 2 > n/2 shows that there exists a 
function C from [0, +oo[ to [0,+oo[ such that, for all e E]0, 1] and all t S 

11^)11^-! ^ ||f/ £ (0)||^_ 1 exp(tC7(M £ (t))), 

where 

M £ (t) := \\U £ \\ LoO ^ 0jt ]. H a-l(^n^ + \\(a, </>)||z / oo([ 0it ].^ CT +2( R n)) • 

It remains to estimate the initial data. By definition, we have 

\\U e (0)\\ Ha -i < 11^(0)11^ + 11^(0)11^. 
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The second term is uniformly bounded by assumption. To estimate the first 
term, recall that 



\z\ 2 -\z'\ 2 



Qa( z i z ')\( z>z >) = ( a e >a y 



where 



Z , \Z' 



l\2\ 



The function Q a is not smooth, but homogeneous of degree a — 1. To 
estimate q 6 at time t = 0, we use the usual product rule in Sobolev space 
and Lemma 077] (applied with F(y 1 , . . . ,y 4 ) = Q CT (yi + %2, 2/3 + «/4)): if 
a" 3, with m = <T — 1 ^ 2, we obtain 



: (0)||^-i 



< 



< 



1 (|« £ (0)| 2 - |a(0)| 2 )]|^ CT i \\QAa £ (0),a(0))\\ Ha 

a e (0),a(0W H ^- 



H <T-l 



The assumption Oq — ao = 0{e) in H s for all s > then implies 
(4.15) sup 11^(0)11^-1 < +oo, 

£610,1] 



hence 



sup \\U £ (0)\\ Ha 
eelo.i] 



< +oo. 



Consequently, since \\u £ e' 



-i<f>/e 



\H° 



Ijjcr-i, the same conti- 
nuity argument as in Section [3] completes the proof of Theorem 14.61 in the 
case <t^3. 

4.4. The case a = 2. For a = 2, we have rei = cr — 1 > re/2 only when 
re = 1. The last point in Lemma 14.71 shows that 



sup \\q 

ee]0,i] 



< +oo. 



We can then proceed as in the case a ^ 3, to prove the second case in 
Theorem 14.61 

Finally, when a = 2 and 2 ^ re ^ 3, recall that we already know that 
for fixed e E]0, 1], u e is global in time, u £ S C(K, H l ). For re = 2, this is 
so since every defocusing, homogeneous nonlinearity is H 1 sub-critical. For 
re = 3, the nonlinearity is H 1 critical, and this property follows from [19J. 
The proof of the estimate is based on an interesting feature of the equation 
for f3 £ (see Proposition I4.3[) . which does not appear in Corollary 14.51 In the 
introduction, we claimed that the previous nonlinear symmetrization of the 
equations implies a local version of the modulated energy estimate. To see 
this, introduce 



,£|2 



+ IV> £ | 2 + kT eC7 i ([0,r e [xI 
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It satisfies an equation of the form dte £ +div(r] £ )+b £ = 0(e £ ), where J b £ = 0. 
Indeed, directly from (|4.10j) . we compute 

d t e £ + div(ve £ ) + 2div(lm(g £ q £ a £ ip £ )) + elm (a £ Aa £ + ^Aip e ) 

= -a\q £ \ 2 divv -Re ((2i(; £ ■ Vv + a 6 V div v)ljF) . 

Hence we have obtained an evolution equation for a modulated energy, which 
yields the desired modulated energy estimate. Gronwall lemma yields 

\\e £ (t)\\ L i {Rn) < l|e £ (0)|| L i (R n ) exp(Ci) . 

Finally, (e e (0)) e is bounded in L 1 (M"). This is obvious for the first two 
terms of e £ . For q £ , a rough estimate yields: 

lk £ (0)|| i2 < je- 1 (V(0)| 2 ~ K0)| 2 )|| i2 II & (« £ (0),a(0))|b 
and the assumption a £ Q — a® = 0(e) in H°° shows that 

sup ||e £ (0)|| L i (R ») < oo. 

0<e<l 



This completes the proof of Theorem 14.61 

4.5. Convergence of position and current densities. As we have 
already mentioned, Theorem 11.31 is a consequence of the first part of (|4.13p . 
Proposition 11.71 follows from both informations in (|4.13|) . Indeed, (|4.13p 
implies the "usual" modulated energy estimate, as in [6l 1391 129] (see also 
[2]). The boundedness of q £ in C([0, T]; L 2 ), and the convexity argument 
IjO]! . yield 

{\a £ (t,x)\ 2 - \a(t,x)\ 2 ) 2 (\a £ (t,x)\ 2 *- 2 + \a(t, x)\ 2 ^ 2 f dx < e 2 . 



sup 

te[o,T] Jm? 
Therefore, 

(4.16) sup f (\a £ (t,x)\ 2 -\a(t,x)\ 2 ) a+1 dx<e 2 . 

te[o,T] JR n 

This yields the first part of Proposition 11.71 along with a bound on the rate 
of convergence as e — » 0. For the current density, write 

Im (eu £ Vu £ ) = \a £ \ 2 V(p + Im (ea £ Va £ ) . 

Since V</> G L°°([0,T] x 1"), ([TO]) yields 

\a £ \ 2 V<t> — >|a| 2 V0 in C([0,r];L a+1 ). 

s— >0 

On the other hand, since a £ is bounded in C([0, T]; H 1 ), we have: 

Im (ea £ Va £ ) =0{e) in C([0, T\;L l ). 

This completes the proof of Proposition 11.71 

Remark 4.10. Since we have used ()4.13j) with k = 1 only, we could also 
refine the statements of Proposition [TT7] when k ^ 2 is allowed in (|4.13p . 
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5. Proof of Theorem 11.91 

To prove Theorem II .91 resume the approach of E. Grenier [26J. His idea 
was to seek 

u £ (t,x)=a e (t,x)e i<t)£< - t ' x)/£ , 

where the pair U £ = (a £ , V^ £ ) is given by a system of the form (|2.4p (with 
E = 0). The point is that the form (|2.4p for this U £ meets all the require- 
ments that we have listed, if and only if the nonlinearity is defocusing, and 
cubic at the origin. In the case of the homogeneous nonlinearity of (jl.ip . 
the only admissible case is then a = 1. The second step of the analysis in 
|26j consists in showing that under suitable assumptions, a £ and 4> £ have an 
asymptotic expansion of the form 

a £ ~ a + ea^ + eV 2 ) + . . . ; (f> £ ~ + e^ 1 ) + e 2 ^ + . . . 

The pair (a,</>) solves (the analogue of) (|1.2I) , Note that because the phase 
£ is divided by e, we need to take into account in order to have a 
point-wise description of u £ : 



u £ ~ ae^ (1) e i0/e . 



Therefore, the rapidly oscillatory phase for u e is given by <j>, and its ampli- 
tude at leading order is given by ae l( ^ m (which does not depend on e). If u £ 
solves 

where Oq satisfies Assumption 11.81 then is given by the system 



ied t u £ + £ -Au £ = f(\u £ \ 2 ) U £ ; uf t=Q = a £ e^ £ , 



dt^ + • V^ 1 ) + 2 Re (aaP>) f (\a\ 2 ) = 0, 



d t a W + Vc/> • Va^ + V0 (1) • Va + -a (1) + -aAc/> (1) = ^Aa, 

^ (1) | t=0 = ; aW| t=0 = a!. 

This coupling shows that <^>W i s a (nonlinear) function of a, <^>, and ai, the 
term of order e in the expansion of the initial data a§. In our case, /(y) = y 
we introduce the system 



£7 . 



(5-1) 



+ V0 • Vc/> (1) + 2aRe (aa w ) \a\ 2a - 2 = 0, 



9ta (1) + V0 • Va (1) + Vc/> (1) • Va + ^a (1) A</> + ^aA^ = % -Aa, 



W| =0 • a«| 
li=0 u ' " \t=0 



ax. 



Lemma 5.1. Let n ^ 1, and let Assumption \1.8\ be satisfied. Then (I5.ip 
has a unique solution ((ft^^a^) in C([0, T*[; H°°(M. n )), where T* is given 
by Lemma \l.2\ . 
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Proof. Again, at the zeroes of a, (|5.ip ceases to be hyperbolic, and 
we cannot solve the Cauchy problem by a standard argument. The strat- 
egy of the proof is to transform the equations so as to obtain an auxiliary 
hyperbolic system for (Vc/^ 1 ), ^4i) for some good unknown A±, depending 
linearly upon The definition of A\ depends on the parity of a. This 
allows to determine a function and next to define a function aW by 
solving the second equation in (|5.ip . We conclude the proof by checking 
that (<p^\ a^) does solve (|5.ip . The first change of unknown consists in 
considering v\ := V^ 1 **. The first equation in (|5.ip yields: 

dm + v ■ Vui + 2cjV Re (\a\ 2a ~ 2 aa {l) \ = -v x ■ Vu, 



where we have denoted v = V0. 

First case: a ^ 2 is even. Consider the new unknown 



Ax := |ar^ 2 Re (aa (1) 

We check that, if (0W,aW) solves dSH]), then 

r a t «i + v ■ V«i + 2cr|a| <T VAi = -«! • Vw - 2cj^iV (|a| CT ) 



(5.2) 



«9 4 yli + v ■ VAi + -|a| CT divui = --V (|a| CT ) • ui - -A x divw 
2 a 2 

+ ^Re (\a\ a ~ 2 aAa) . 



This linear system is hyperbolic symmetric, and its coefficients are smooth 
since a G 2N and a,v G C°° ( [0, T* [; iif °° (IR n ) ) , from Lemma E2 In partic- 
ular, uniqueness for (|5.ip follows from the uniqueness for (|5.2p . Note that, 
since a - 2 G 2N, 



KA)Lo = (0,|a r- 2 Re(a ai)) G " " 2 

Therefore, ([521) possesses a unique solution in C°°{[0,T*[; H°°(W 1 )). We 
next define G C°°([0,T*[: -£P°(IR n )) by 

0( 1 )(t,j C ) = - / (w(t,x) • wi(r,x) + 2<T|a(T,x)| <T ^i(r,x)) dr. 



Then (V^ 1 ) — v\) = 0, therefore v\ = and hence (j)^' satisfies 

d t (/) {1) + u • V0 (1) + 2<r|a|Mi = 0, | i=Q = 0. 

Once this is granted, we can define aW G C°°([0,T*[: H°°(R n )) as the 
unique solution of the linear equation 

<9 t a (1) + v ■ Va (1) + V0 (1) • Va + ^a (1) divv + ^aA^W = ^Aa, 
a«| f=Q = ai. 
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By construction, A\ and \a\ a 2 Re(aa( 1 )) solve the same linear equation, 
where (jP^ is viewed as a smooth coefficient. Therefore, these two functions 
coincide, and (<p^\a^) solves f)5. 1[) . 

Second case: a is odd. In this case, a = 2m + 1, for some m € N. We 
consider the new unknown 



\a\ 2m aV. 



We check that (v\ , Ai ) must solve 
' d t vi + v ■ Vui + 2cj Re (|a| 2m aV^i) = -«i • Vtj - 2a Re (AiV {\a\ 2m a)) 



d t A 1 + v • VAx + -|a| 2m adivwi 



-— divu 
2 



|a| 2m Va- 



+ 2 |a 



2m 



Aa. 



□ 



We can then conclude as in the first case, by considering (vi, A\, A\) 
Theorem 11.91 follows from: 

Proposition 5.2. Let n < 3, and Zei Assumption \1.8\ be satisfied. Set 
a := ae l W . Then for any T G]0, T*[, there exists e(T) > such that 
a £ G C7([0,T];F°°) /ore G]0,e(T)], and 

||a £ — a|| ioo([o,T];H fc ) = ^( e )> 
where k is as in Theorem \J7t 



Proof. Since the proof follows the same lines as the proof of Theo- 
rem [L3l we shall indicate its main steps only. Denote 

r £ = a £ -a ; a« = a« . 
From (jl.2p . (jl.6p and (|5.ip . we see that r e solves 

dtr £ + v ■ Vr £ + -r e divv — i-Ar £ = i-Aa — iS £ , 
2 2 2 

rj* t=0 = Oq - a = ea\ + O (e 2 ) , 
where the term S e is given by: 



- | \<r 



2<T 



2o-a\a\ 2a - 2 Re(aa {1) 



We check that for all s ^ 0, we have, in fP(R n ): 



2cr 



a + ea 



(i) 



2cr 



a e + 2 CT r e |5| 2,J - 2 Re fe 1 )) +0(e). 



The last term should be viewed as a small source term. The second one is 
linear in r £ , and is suitable in view of an application of the Gronwall Lemma. 
There remains to handle the first term. At this stage, we can mimic the 
approach detailed in Introduce the nonlinear change of unknown: 



1 



-B n 



a + ea 



(i) 



g £ = G 



a + ea 



(i) 
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where B a and G a are denned in Notation 14.11 We check that (r £ , Vr e , ) 
solves a system of the form (|4.1ip . plus some extra source terms of order 
0(e) in H s (W n ). We also note that the initial data are of order 0(e), from 
Assumption 11.81 

\\(r £ ,Vr% =Q \\ Ha =0(e), Vs ^ 0. 

We also have 

ll5l t =olU-i = °( £ )> 

where k is as Theorem 14.61 

Following the approach of £JH the proposition stems from Gronwall 
lemma. Note also that the time T can be taken arbitrarily close to T*, 
by the usual continuity argument, since we now have an error estimate that 
goes to zero with e. □ 

To conclude this paragraph, we note that unless ao is real valued and 
ai G iffi, one must not expect o = a. Indeed, we see that 

^0 = ; d t ^l Q = -2(jRe(aoa 1 )|a | 2,J - 2 . 

So in general, </>W ^ 0, and a ^ a. On the other hand if ao is real-valued, 
then so is a. In this case, 

Im (oAo) = 0, 

and Re(aoW)) solves an homogeneous linear system. Therefore, if 

Re(aa (1) ) = at time t = 0, then c^^ 1 ) = 0. 

6. Further remarks 

The following remarks serve to clarify some features of the systems we 
produced. 

6.1. Regularity of the initial data. It is a matter of routine to 
extend the previous analysis to the case where the initial data belong to 
H s (W l ) with s < +oo large enough. 

6.2. Introducing an external potential. To treat a possibly more 
physically relevant case, one might want to consider (jl.l|) with an extra 
external potential: 

ied t u £ + — Au £ = Vu £ + \u £ \ 2a u £ ; u £ t=0 = age^ o/e , 

where V = V(t, x) is real-valued, and possibly time-dependent. As noticed 
in |11| , it is sensible to consider an external potential V and an initial phase 
4>o which are smooth and sub-quadratic: 

d%V G C(M; L°°(]R ri )), d a <f> G L°°(M. n ), Va G N n , \a\ ^ 2. 

This includes the case of the harmonic oscillator, commonly used in the 
theory of Bose-Einstein condensation ( |27| ). The main remark in |llj is 
that the introduction of this assumption does not deeply change the analysis. 
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Indeed, we can resume the analysis of (ll.lj) : introduce the solution to the 
standard eikonal equation 

1. l2 I 

dt&eik + ^x^eikl +V = ; 4>eik\ t=Q = 4>Q ■ 

Decomposing the phase <p of the above quasi-linear analysis as 

<t> = 0cik + <P, 

and seeking (ft in Sobolev spaces, we see that the extra terms appearing 
after this sort of linearization can be treated like semi-linear terms. There- 
fore, mimicking the above computations, and using only extra perturbative 
arguments, it is easy to adapt Theorems 11.31 and 11.91 to this case. 

6.3. About conservation laws. Recall some important evolution laws 
for ([LTD : 

Mass: ^\\u £ (t)\\ L 2 =0. 

Energy: | Q| k V^||| 2 + -^||^||-t 2 2 ) =0. 
d f 

Momentum: — Im / u £ (t,x)e'V x u £ (t,x)dx = . 
at J 

d (\ t 2 \ 

Pseudo-conformal law: — ( -II J £ (t)u £ \\ ? 2 H l|w £ ||?2i"+2 

dt \2 ^ a + 1 J 



where J e (t) = x + ietV x . These evolutions are deduced from the usual 
ones (e = 1, see e.g. [13|, 133] ) via the scaling tp(t,x) = u(et,ex). Writing 
u £ = a £ e l( ^l £ , and passing to the limit formally in the above formulae yields: 

d 



dt Ht)h*=o. 

j t J (^\ a (t,x)\ 2 \V ( l ) (t,x)\ 2 + -^\a(t,x)\ 2 ^dx = 0. 

d f 

— / \a(t,x)\ 2 V(f)(t,x)dx = 0. 



dt 

| J (^\( x -tV^t,x))a(t,x)\ 2 + -^\a(t,x)\ 2 ° +2 )dx 



a + 1 

Note that we also have the conservation (|12j): 
d 



' (2-na) I \a(t,x)\ 2a+2 dx. 



^Re / v £ (t,x)J £ (t)u £ (t,x)dx = 0. 
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which yields: 

^- I (x - tV<j)(t,x)) \a(t,x)\ 2 dx = 0. 



dt . 

All these expressions involve only (|a| 2 , V</>) = (|a| 2 , V<^>). Recall that if we 
set (p,v) = (|a| 2 , V0), then (|1.2|) implies 

^ f + « • + V = ; v lt=0 = V0 O , 

| d t p + div (pu) = ; p\ t=Q = |a | 2 . 

Rewriting the above evolution laws, we get: 

d 

dt J J^n 



p(t, x)dx 



^ y p(t, x)v(t, x)dx 



dt . 

1 / Ql^-Mi^))| 2 p(*^) + ^;P(i^) a+1 V^ = 

(6.2) = -^-y(2-ncr) f p(t,x) a+1 dx. 

d f 

— / (x — tuft, a;)) p(t, x)dx = . 
dt J 

We thus retrieve formally some evolution laws for the compressible Euler 
equation (|6.ip (see e.g. |32L I37| ). with the pressure law p(p) = cp a+1 . 

6.4. About global in time results. We point out that the solution 
to (|1.2p must not be expected to be smooth for all time: the time T* in 
Lemma 12.21 is finite in general. Recall that (p,v) = (|a| 2 ,V0) solves (|6,ip . 
Theorem 3 in |30| (see also |37j ) implies that, if V^o and |ao| 2 are compactly 
supported, then the life span T* in Lemma [2. 2 1 is necessarily finite. Note that 
these initial data can be chosen arbitrarily small: the phenomenon remains. 

Proposition 6.1. Let n ^ 1 and a ^ 1. For all initial data (ao,0o) £ 
C 2 (lR n ) with compact support, there does not exist (a,(f>) £ C 2 ([0, +oo[xM n ) 
satisfying the Cauchy problem (|1.2p . 



A word of caution: because of one technical assumption in the defi- 
nition of regular solution in [30] . Theorem 3 in |30j does not apply di- 
rectly. Yet, one can prove our claim by combining the proof of Lemma 12.21 
with the approach in |30| . Indeed, recall that U := (a°",V0) satisfies 
dtU + Aj(U)djU = where the Aj's are n x n matrices linear in their 
argument. Therefore, the proof of Theorem 2 in [30 j shows that U is com- 
pactly supported, and so is (p, v) := (|a| 2 , V0), with support included in the 
support of (|ao| 2 , V<^q). And this is the only point which requires the above 
mentioned technical assumption. 



GEOMETRIC OPTICS FOR NLS 



27 



Note also that the proof of this result in [37] relies on the evolution law 
for the total pressure 



This approach is very similar to the Zakharov-Glassey method [381 125] , 
which yields a sufficient condition for the finite time blow-up of solutions to 
the focusing nonlinear Schrodinger equation. As noticed by M. Weinstein 
[36] . the identity used by Zakharov, and generalized by Glassey, follows 
from the pseudo-conformal law, along with the conservation of energy. For 
a ^ 2/n and a defocusing nonlinearity, this approach yields an upper bound 
for the L 2 -norm of xu, the momentum of u. When this upper bound may 
become negative, finite time blow-up occurs. 

In the present context, the nonlinearity is defocusing, but the idea is 
similar. Note that (the generalized version of) (16. 2j) is the key ingredient 
in the proof of Z. Xin [37J (Z. Xin considers Navier-Stokes equations). Ex- 
panding (16. 2|) . and using the conservation of energy, we recover an upper 
bound for (|6.3p which goes to zero as t — > oo. But so long as v remains 
bounded, (16. ip is an ordinary differential equations for p, thus contradicting 
the upper bound for (|6.3j) . unless v ceases to be smooth in finite time (see 
[37] for the details). 

6.5. About focusing nonlinearities. The main feature of the limit 
system we used is that it enters, up to a change of unknowns, into the 
framework of quasi-linear hyperbolic systems. This comes from the fact that 
we consider the defocusing case. Had we worked instead with the focusing 
case, where +|u| 2<J u is replaced with — |u| 2cr u, the corresponding limit system 
would have been ill-posed. We refer to [31] . in which G. Metivier establishes 
Hadamard's instabilities for non-hyperbolic nonlinear equations. 

As an example, consider the Cauchy problem 



The following result follows from Hadamard's argument (see |31j). 

PROPOSITION 6.2. Suppose that {4>,a) in C 2 ([0, T] x E) solves (|6.4j) . If 
<f>o(x) is real analytic near x and if ao(x) > 0, then ao(x) is real analytic 
near x. Consequently, there are smooth initial data for which the Cauchy 
problem has no solution. 

This shows that to study the semi-classical limit for the focusing ana- 
logue of (jl.ip , working with analytic data, as in [221 135] , is not only conve- 
nient: it is necessary. 



(6.3) 




(6.4) 
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